Cooperative Electronic and Phononic Mechanism of the High Temperature 

Superconductivity in Cuprates 
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In conventional superconductors, phonons glue two electrons with opposite spins to form Cooper 
pairs and condensation of these pairs leads to the superconductivity. Identifying the underlying 
mechanism of the high temperature superconductivity in cuprates is among the most important 
problems in physics. Even quarter of a century after the first report of high temperature supercon- 
ductor by Bednorz and MuUer in 1986, there is still no general consensus on the pairing mechanism 
of superconductivity in these materials. So far, many theories have been developed to explain the 
exotic properties of cuprates, but they can explain only a limited number of experiments. In this 
article, we present a new pairing mechanism that incorporates both strong correlation and phonon 
mediated interaction on an equal footing to produce superconductivity. In this framework, strong 
correlation and anti-ferromagnetic interaction between electrons, create RVB pairs and phonons 
provide the phase coherence between these RVB pairs. Both of these are required in this approach 
to obtain the superconductivity. This approach resolves three limitations of the U(l) slave boson 
method. We achieve a better estimation of Tc, we only predict vortices and the linear T coef- 
ficient of the superfiuid is not sensitive to the doping. This formalism provides a framework that 
connects Anderson's idea of preformed Cooper pairs and phonon based theories. 



I. INTRODUCTION 

One scenario for the high temperature superconduc- 
tivity in cuprates [l| is the preformed Cooper pairs idea 
that was first proposed by Anderson 2]. This theory is 
based on the strong correlation effects and it does not 
incorporate phonons in the pairing mechanism. On the 
other hand, many experiments have been reported indi- 
cating the importance of electron phonon interaction in 
understanding the physics of high Tc cuprates, such as 
the oxygen isotope effect on both the transition temper- 
ature of superconductivity and the London penetration 
depth 

Observation of the strong isotope effect in cuprates has 
led many authors to employ the strong limit of the elec- 
tron phonon interaction as the primary cause of the su- 
perconductivity in these materials [9| . For example, some 
workers have applied bipolaron theory of superconduc- 
tivity to the high temperature superconductivity prob- 
lem. This theory requires the breakdown of the Migdal- 
Eliashberg theory and is based on the non-adiabatic limit 
of the electron phonon interaction, where phonons have 
a much larger energy than electrons. Experimentally, a 
typical energy of electrons is around the exchange energy 
J ~ 130 meV 11(1 . while the energy of optical phonons is 
40 — 70 meV fll[. On the other hand, the breakdown of 
the Migdal-Eliashberg theory 12| in any phonon based 
theory is crucial, because isotope experiments in cuprates 
are very different from conventional superconductors that 
are explained by the BCS theory 13| and its generaliza- 
tion, i.e. Migdal-Eliashberg theory. 

Some authors emphasize on the importance of the 
microscopic inhomogeneity, charge and spin stripes in 
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understanding the mechanism of the superconductivity 
[T3 . [l5| . The idea is that the phase segregation can save 
kinetic energy of holes and the exchange energy of spins. 
Stripe models are based on the competition between this 
phase and other states of matter. It is noteworthy that 
static stripes have been observed only in La2SCu04 fam- 
ily near x = i, not in YBCO family, but dynamical 
stripes are expected even away from x — ^ [16]. The 
idea of Anderson can be implemented using the slave 
particle method [10, 17-20]. The U(l) version of this 
method explains many basic properties of cuprates suc- 
cessfully. However, this method has several limitations. 
For example, the U(l) slave boson method overestimates 
the transition temperature of the superconducting state 
by an order of magnitude. Although in experiments only 

vortices have been observed, this method predicts an- 
other — vortices as well. For clean d-wave superconduc- 
tors, the superfiuid density decreases linearly with tem- 
perature up to the leading order. Experimentally the 
coefficient of this linear term is not very sensitive to the 
doping, while the U(l) slave boson treatment predicts a 

doping dependence, where x is the doping percentage. 
In this paper we argue that phonons can mediate attrac- 
tive interaction between holons and therefore we find a 
paired holon state. We show that this assumption re- 
solves the three mentioned limitations of the U(l) slave 
boson approach. 

Here we extend the Anderson theory of high Tc which 
is based on the resonating valence bond(RVB) state, by 
adding Holstein Hamiltonian such that the model Hamil- 
tonian include the effect of electron phonon interaction. 
RVB state is the superposition of all possible singlet 
states between any two sites. The idea of RVB state can 
be quantified using the slave boson method, and natu- 
rally leads to the spin charge separation in two spatial 
dimensions. The low energy excitations of this state are 
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described by the charged spinless quasiparticles which are 
called holon and are treated as bosons, and the neutral 
spin 1/2 quasiparticles which are called spinon and are 
treated as fermions. This can be implemented by writ- 
ing the creation operator of physical electrons cj ^ as the 
product of a fermionic operator fj^ and a bosonic oper- 
ator hi. Therefore ^ ^ fla^i- This assumption comes 
along with the nondoubly occupancy constraint due to 
very large onsite repulsion between electrons in cuprates. 
In terms of slave particles fj^fi.t + fi^.fi,i + h\hi = 1 at 
every site. Therefore at each site, we have either a holon 
or a spinon. Using this relation, the density of holon gas 
is set by doping x, and the density of spinous is equal 
to the density of electrons 1 — a;. Spinon pairing forms 
Cooper pairs in the system and holon condensation pro- 
vides the phase coherence for these pairs and the super- 
fluid density of the superconducting state is controlled by 
the condensation fraction of holons. Therefore, to obtain 
superconducting phase, both holon gas and spinon gas 
should condense. In case of spinous, we need a pairing 
mechanism due to their fermionic nature. Strong antifer- 
romagnetic interactions in cuprates can result in such a 
pairing potential and as a result the transition tempera- 
ture of spinon paired state is controlled by the exchange 
energy J which is quite large. This large gap in excita- 
tion spectrum of spinous is the origin of the pseudogap 
phenomenon in cuprates. 

On the other hand, due to the bosonic nature of holons, 
they do not need a pairing mechanism and they can un- 
dergo Bose Einstein condensation(BEC). Although this 
can be in general true, we discuss in the following sec- 
tions that the phenomenology of cuprates supports the 
holon pair condensation scenario. This assumption im- 
mediately resolves several limitations of the Anderson 
theory. 

In order to pair holons, like in the BCS theory, we need 
a pairing mechanism. In our theory, electron phonon in- 
teraction can cause such a pairing instability. As we have 
shown in Ref. 0, this assumption can also explain the 
observed phonon experiments successfully. In summary 
in this framework spinous form RVB pairs and holons are 
also paired by phonons. Both of these are required for 
the superconductivity in cuprates. 

This approach, gives a phase diagram similar to the 
phase diagram of the conventional U(l) slave boson treat- 
ment of the t-J model, except that the superconducting 
state has Z2 gauge symmetry, and we obtain a better es- 
timation of Tc, closer to the ex per imental data(see Fig. 
[T|) . We have also shown in Ref. [20| , that at least at half 
filling, the antiferromagnetic order in the Hubbard model 
can be achieved by including a nonzero triplet component 
in the pairing amplitude besides the singlet component 
[20| . Now let us briefly comment on the experimental 
consequences of our theory. In appendix, we present more 
details and derive equations. 




FIG. 1: Phase diagram. — Schematic phase diagram of 
cuprates from our theory. Xopt is the optimal doping percent- 
age which is typically around %15. It corresponds to the op- 
timum transition temperature which is typically around 100 
Kelvin. SC is the superconducting state in which we obtain 
spinon pair condensation and holon pair condensation. PG 
is the pseudogap phase where spinous are paired but holon 
gas in not condensed. FL is the Fermi liquid phase where we 
obtain single holon condensation but no spinon pair conden- 
sation. SM is the strange metal in which neither holons nor 
spinous condense. In this approach, spinous form RVB pairs 
due to anti-ferromagnetic interactions and holons are paired 
by phonons. Both of these are required to obtain the super- 
conducting phase in cuprates. We have not shown AF order 
in this phase diagram, however we have shown in Ref. [l^] 
that it can be achieved by adding a nonzero triplet component 
to the pairing amplitude near half filling. 



II. BETTER ESTIMATION OF Tc 

In the underdoped cuprates, both the slave boson the- 
ory and the phase fluctuation studies of superconduc- 
tors predict that the transition temperature of the su- 
perconducting state is controlled by the superfluid den- 
sity. Using loffe-Larkin recombination formula (see Ap- 
pendix C for derivations and discussions) , it can be shown 
that at small dopings, the superfluid density is mostly 
determined by that of holons. Therefore Tc is deter- 
mined by the BEC transition temperature T^^^ , of holon 
gas. In the single holon condensation theories we ob- 
tain T^^^ — Although this gives a good doping 

h 

dependence, it is an order of magnitude higher than the 
superconducting transition temperature [2l|. In the pair 
condensation scenario, we obtain the following expression 
for Tc (see Appendix B for derivation) 

/ ~\ T^ 
TP (x,V]^ ^ (1) 

^^"^ ^ -ln(2yiV(0)) ^ ^ 
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where V = VA^, where tAg is the pseudogap energy, is 
the renormahzed couphng constant of the phonon me- 
diated holon holon attraction, and A^(0) = ^ is the 
density of states. In Ref. [1], we had to chose V = 12^, 
where is the energy of optical phonons, to fit theory 
and experiment. In cuprates, 12wg ~ t and therefore 
V ^ A^ (Tc) ^ 1. Because of the extra factor in the de- 
nominator, T^^^ is much smaller than TJ^^ and is closer 
to the experimental data. 



III. VORTICES 

In the single holon condensation scenario, two kinds of 
vortices are allowed, and — . Studies show that the 

' 2ec ec 

energy of the latter is much smaller, so it should be more 
stable and visible in experiments. In experiments how- 
ever, only 2^ vortices have been observed [l^. Absence 
of ^ vortices, challenges the assumption of single holon 
condensation in cuprates. In the holon pair condensa- 
tion scenario we always obtain vortices from both 
(hh) or (/f/j,) order parameters. This in fact reflects the 
Z2 structure of our theory. 



tides have effective charge — e (see Appendix D for de- 
tails) . So we obtain a = 1 in this case, in consistent with 
the linear temperature coefficient of superfluid density 
measurements. 

In Ref. Q, we have investigated the isotope effect on 
the superfluid density and the transition temperature of 
the superconducting phase by considering the effect of 
electron phonon interaction and have achieved a good 
agreement between theory and experiment. 

It is worth mentioning that boson gas in a purely at- 
tractive potential is unstable. It collapses and phase sep- 
aration happens in that case. However in our model, we 
deal with holons which are hardcore bosons, i.e. there is 
an infinite on-site repulsion between them. The phonon 
mediated attraction is also screened by the presence of 
spinous. These two can make the paired state of holons 
stable. We speculate that under some certain conditions, 
holons may become meta-stable and form stripes due to 
the phonon mediated attraction. Moreover as we men- 
tioned earlier, the phase segregation saves the kinetic 
energy of holes and exchange energy of spinous. This 
can further stabilize stripe order phase. Therefore, the 
phonon mediated attraction between holons may enhance 
the stripe formation. 



IV. LINEAR T COEFFICIENT OF THE 
SUPERFLUID DENSITY 

Another important limitation of the single holon con- 
densation scenario, is the calculation of the linear tem- 
perature coefficient of the superfluid density. This sce- 
nario predicts a parabolic doping dependence behavior, 
while experimentally it has a weak dependence on the 
doping percentage [2l], [l^. The reason is that within 
this assumption, the current carried by quasiparticles 
is j = aeVp , where is the Fermi velocity of nodal 
quasiparticles and a ~ a; (see Appendix D for details). 
Since holons are charged particles, they couple to both 
the external gauge field (Aext)i and the induced inter- 
nal gauge field Aint. Within single holon condensation 
scenario and using loffe-Larkin formula, it can be shown 
that when both spinons and holons condense, we have 
Aint ~ —xAext- Since spinons are electrically neutral, 
they only couple to the internal gauge field and therefore 
they see —xAext, so their effective electric charge is —xe. 
Now we can estimate the value of a by computing the 
Green's function of real electrons. In the single holon 
condensation scenario, the diagonal part of the Green's 
does not respond to the gauge field, because it depends 
on (h^) {h). Therefore only spinons couple to internal 
gauge field which is by factor of —a; smaller and this is 
why a — X in this case. In the pair condensation scheme, 
since (h) = 0, the diagonal part of holons Green's func- 
tion depends on (^h^^hi^. So it responds to the external 

as well as the internal gauge field. From convolution it 
can be checked that the real electrons respond to the 
whole external electromagnetic field and there quasipar- 



V. CONCLUSION 

In this paper, we have extended the Anderson theory of 
high Tc to take phonons into consideration. Phonons are 
engaged in the pairing mechanism by mediating attrac- 
tive interaction between charged spinless quasiparticles 
(holons). This attraction destabilizes the single holon 
condensation state and gives rise to the paired holon 
state. Assuming the paired holon phase immediately re- 
solves several limitations of the U(l) slave treatment of 
the Anderson theory. First of all, we achieve a better es- 
timation of the superconducting transition temperature. 
Moreover, since both spinons and holons are paired, we 
only find vortices. Finally we showed that the effec- 
tive electric charge of nodal quasiparticles is —xe when 
we have single holon condensation, so they carry xeVp 
current, while the the effective electric charge of nodal 
quasiparticles is — e and they carry ew^ when holons con- 
dense in pairs. We have shown that in the latter case, 
the linear temperature coefficient of the superfluid den- 
sity is almost independent of doping, consistent with ex- 
perimental observations. 
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APPENDICES 

In the following, we provide detailed derivations of for- 
mulae in the main text concerning the expression for Tc, 
loffe-Larkin formula and Linear temperature coefficient 
of superfluid density. In Appendix A, we provide a brief 
introduction to the slave boson method. In Appendix B, 
we derive an expression for the transition temperature 
of the superconducting state. In Appendix C, we derive 
loffe Larkin formula which serves an a powerful tool in 
relating physical quantities to the corresponding proper- 
ties of the slave particles. In Appendix D, we present an 
argument to calculate the linear temperature coefficient 
of the superfluid density. 

APPENDIX A: METHOD 

Let us start from the t-J model as our starting point. 
It is believed that this Hamiltonian captures the essen- 
tial physics of the strongly correlated systems. We finally 
add the Holstein Hamiltonian to take the effect of elec- 
tron phonon interaction into consideration. t-J model is 
defined as the following 

where Pq is the Gutzwiller projection operator that 
removes doubly occupied states. Within the slave bo- 
son formalism, electrons can be decomposed as c| g. = 
f} ijhi along with the physical constraint on each site: 

h\hi + fj afi,(y — 1 which implements the Gutzwiller 
projection. / particles are fermions and we call them 
spinon and h particles are bosons and we call them holon. 
Spinon corresponds to a state with only one electron and 
holon to an empty site. The definition of the projected 
electron operator along with the constraint, says that we 
have always one slave particle per site and two slave par- 
ticles cannot sit on the same site. On the other hand, 
whenever cl ^ acts on an empty site, it annihilates one 
holon and creates a spinon with spin a. We cannot act 
further on the resulting state by cj since this operator 
has to kill a holon, but there is no holon anymore at that 
site. If we act Ci^a- on a site that contains a spinon with 
spin cr, the operators annihilates the spinon and creates 
a holon at that site. So by acting projected electron op- 
erator we always annihilate one type of slave particle and 
create another one and therefore the number of slave par- 
ticles at each site is conserved. Now we can rewrite the t- 
J model in terms of the new slave particles. Within mean- 
field approximation and by using Hubbard-Stratonovic 
transformation, we can decouple spinous (spin sector) 



from holons (charge sector) and we obtain the following 
effective Hamiltonians for each sector 

^^/i = - E ^Xsh\hj - '^fifihlhi (3) 

<i,j> i 

- {J 12) A, [jyi - fljl^ + H.cM) 

<i,3> 

where the following notations have been used 

X. = {hl,h,) (5) 
^^-(EfllJ-^) (6) 

^si^,J) = l{fyh-fifl,) (7) 

At low temperatures, most of holons occupy the 
groundstate with momentum k = 0, therefore Xh ~ x. 
This model has been extensively studied in the literature 
and it is well known that this model leads to the d-wave 
pairing symmetry for spinous [lo| . i.e. (±x) — 
and A,. (±y) = - A,. 

Now let us consider the electron-phonon interaction. 
Since the typical energy of electrons is around J and is 
much larger than , we can apply the BCS theory in our 
case. Within BCS theory, electron phonon interaction, 
leads to the following pairing term: 

- E ^fe.fe' < 4',tC-/c'4 > c-fe.;cfe,t (8) 
fc.fc' 

By translating the above term to the slave boson lan- 
guage in real space, and using the mean-field approxi- 
mation, we can substitute Ci_iCj^f = fi,ifj,-thlhj by the 
following terms: 

Ah (i, j) /u/m + As (i, j) hlh] - A, Ah (i, j) (9) 

where A^ (hj) ~ (hihj) ^ x (doping). Now assuming 
a very short range interaction, we obtain the following 
effective interaction: 

H[_, = -VAl ^ A, {^,J).fyl + B.C.. (10) 

<id> 

HU, = -VAl J2 (^'^■) f'l'^] + H e (11) 

<i,3> 

2 

in which V = ^^'^ 7o is the bare electron phonon 

E 

interaction coupling constant, and lOj^ is the energy of 
optical phonons. Let us assume that VAf^ ~ Vx"^ <C J, 
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so we can neglect this phonon mediated pairing term 
and therefore, the d-Wave nature of the spinons does 
not change. From the above we see that the couphng 
constant of phonon mediated spinon spinon attraction 
is renormahzed by factor and that of holons by 
due to strong correlation effects. It is easy to show that, 
y cx 7^, where 7 is the electron phonon coupling con- 
stant. Therefore, we can interpret these renormaliza- 
tion factors as the renormalization of the coupling con- 
stant of spinon-phonon interaction to Afiji^ q and that 
of holon-phonon interaction to Agjj^ q. In Appendix C, 
using loffe-Larkin recombination formula, we present an- 
other physical argument to justify this result. Now let 
us focus on the charge sector (holons). We have studied 
the effect of holon-phonon interaction on the holon mass 
in Ref. Since we treat charge sector as a Bose gas, 
we have shown the holon-phonon interaction is in the 
non-adiabatic limit and the small polaron picture can be 
applied. The mass of holons enhances by an exponential 
factor and we have: to^ (T) = 2x^(t) ~ ^'^^''^h: where 

9HT) = ^. 



APPENDIX B: CALCULATION OF n 

Taking the mentioned mass renormalization into ac- 
count, the following continuum model describes the low 
energy physics of holons: 

Hb = J2 + ImI) hlhk ~ VAh hlh-k + H.c. (12) 



where efc = A^ = Y.k ^kh^k, and V = VA^ It 

has been shown that in 2D, no matter how small V is, we 
have a bound state of bosons and then the double conden- 
sation is energetically favorable. Using the Bogoliubov 
transformation, we can find the energy eigenvalue which 

is Ek = y {ck + ImI)^ - (vAh^ and the energy eigen- 
vectors. On the other hand, we should choose /z such 
that X = jfj^k^l^k- The two mentioned constraints 
lead to the following self-consistency equations: 



2Ek 



(13) 



(14) 



in which n^^ {Ek,T) = exp(_Efc/T)-i Bose Ein- 

stein distribution function. The first constraint can 
be solved exactly by ■^J2k ^ 1^ ! '^^ ^ ^(0)/ de-k^ 
where N (0) = . After integration we obtain: 



E, = ^^j\n\^-{vA,,\ = -2T\n\^^^±f^ \ where 



y = exp (^ ^^^2T^° ^ and To = f^. At T we al- 
ways have -Eg = in consistent with Hugenholtz-Pines 
theorem. For small values oiV , we have : 



Eg (T) = Texp 



2Ta 



2T 



(15) 



which is diminishing very rapidly and we can neglect 
it up to the first order approximation. Therefore, at 
small enough temperature energy excitations are sound- 
like an are of the form Ek — ck where c — ^ Mr . At 

V '1 

TP, Afc = and therefore Ek = ek + and Eg = 
\n\ = TPexp( l^lyp ) ~ TPexp(-^). The second 
constraint can be written in the following way: 



N{0)V 
4 

N{0)V 
TP = 



coth 



2TP 



■de 



2TP 



(16) 



2ttx 



ln27V(0)y 



-m^ In 



(18) 



Using the paired state and assuming gapless excita- 
tions {Eg = 0), we can show that up to the first or- 
der, (hlfhjfl'^ — Xpc, where Xpc is the ground-state 

macroscopic occupation number (condensation fraction) , 
though (hi^t) — 0. If we replace every operator by minus 
itself the effective Hamiltonian and all order parameters 
remain the same, and therefore we have Z2 gauge free- 
dom and the low energy theory is described by a Z2 gauge 
theory. As long as T < T^, the condensation fraction is 
nonzero and T — 0, Xpc — x, where x is doping (number 
of holons). 



APPENDIX C: lOFFE-LARKIN FORMULA 

One important question that should be addressed is 
how to relate the physical quantities to the correspond- 
ing quantities of spinons and holons? For example, given 
the conductivity of spinons and holons, what is the con- 
ductivity of real electrons? To answer this question we 
should note two things. First of all, since at each site we 
should have one slave particle, therefore if one spinons 
hops from site i to site j, it should be accompanied by a 
hopping of one holons from site j to site i. So we conclude 
that the current carried by spinons is equal but opposite 
to the current carried by holons and they add up to zero. 
Since electron operator at site i with spin a, is written 
in terms of slave particles as = fj ^hi, it is invariant 
under U(l) gauge transformation, provided spinons and 
holons carry the same charge under this transformation. 
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On the other hand, we assume spinons to be electrically 
neutral and assign +e electric charge to holons (we could 
also assume neutral holons and assign — e electric charge 
to spinons). By scaling the internal gauge field we can 
assume Cmt — e. As we discussed we should satisfy the 
following constraint 



Js + A = 



Since only holons carry electric charge we have 



(19) 



(20) 



Above two equations tells us that nonzero external 
gauge field (electromagnetic field) induces internal gauge 
field. Because we have 



Js — eOgEint 

Jh — eOh {Eext + Eint) 



Solving Eq. [19] 



Eint 



- Ep 



and as a result Jph = e E^xt- Equivalently 



Jph 



(Jh(Js 
(Js + (Jh 



(21) 
(22) 



(23) 



(24) 



Now let us consider two important cases. Fermi liquid 
phase and the Superconducting phase. 



A. Fermi Liquid Phase 

In the Fermi liquid (FL) phase, holons are condensed 
but spinons are not. So we have {h) ^ and Ag — 0. 
Since in this case, ah ^ (Js, we have Eint — —E^xt, 
aph — (Js, and Jph — eUsEext- Now let us define the 
effective electric charge of spinons and holons as 



Jph ~ —Gs(JsEext — eh(JhEe. 



(25) 



Therefore in the Fermi liquid we have ~ — e and 
therefore they have nonzero overlap with physical elec- 
trons and carry the same charge. On the other hand 

— and we can safely assume that holons are electri- 
cally neutral in the FL phase. This can bee directly seen 
from ~ {h) f}^. One important result of this argu- 
ment is that in the absence of the pseudogap (i.e. when 
As = 0), condensed holons do not couple to phonons, 
since phonons only couple to electrically charged quasi- 
particles. In other words, phonons create local electro- 
magnetic field and this field induces another (internal) 



gauge field. Holons couple to the sum of these two fields. 
When As = 0, but holons condense {i.e. at low enough 
temperatures) these two fields cancel out each other and 
as a result holons do not couple to phonons. 



B. Superconducting Phase 

In the superconducting (SC) phase , both holons and 
spinons are condensed. So we have A^ = {hh) ^ 
and As ^ 0. In this case, ^ = ^^A^ where pc,ft. and 
Pc,s is the condensation fraction of holon and spinon 
gas respectively. At zero temperature all holons and 
spinons condense and therefore ^ = and we have 

Eint — -xEext, (Jph — xas, and Jph ^ xeasEext- We 
can compute the effective charge of spinons and holons 
and we obtain Cs — —xe and e/i ~ (1 — x) e respectively. 
Therefore in the superconducting state, spinons only re- 
spond to the X fraction of the electromagnetic field. Since 
holons are effectively charged quasiparticles in this state, 
(or since the internal gauge field does not cancel out 
the local electromagnetic field of phonons completely) 
they couple to phonons. The larger pseudogap value, the 
stronger interaction between holons and phonons. This is 
an intuitive way to justify our recipe for renormalization 
of the electron phonon coupling constant from the brae 
value 7, to the effective value As7. Using this expression 
and since As decreases with doping, we can easily explain 
why the isotope effect is a decreasing function of doping 
as well. Moreover, in the overdoped region. As = at 
the transition temperature and therefore holons do not 
interact with phonons and as a result we do not expect 
isotope effect on the Tc. However at T = 0, As 7^ even 
in the overdoped region and the mass of holons enhances 
in an isotope dependent way and this explains the non- 
vanishing isotope effect on the superfluid density and the 
London penetration depth in this region. 



APPENDIX D: LINEAR T COEFFICIENT OF 
THE SUPERFLUID DENSITY 



Another important limitation of the single holon con- 
densation scenario, is the calculation of the linear tem- 
perature coefficient of the superfiuid density. This sce- 
nario predicts a parabolic doping dependence behavior, 
while experimentally it has a weak dependence on the 
doping percentage pT!, '55]. The reason is that within 
this assumption, the current carried by quasiparticles is 
j — aeVp , where is the Fermi velocity of nodal quasi- 
particles and a ^ X. Lee and Wen have shown that the 
linear temperature dependence of a d-Wave superconduc- 
tor is given by the following expression 



Ps {T) 



m 



X 

m 



21n2 



V2 



(26) 
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where V2 is the velocity of the d-wave SC quasiparti- 
cles in the direction perpendicular to v^^ . Now let us give 
a simple argument on how to compute a. Since holons 
are charged particles, they couple to both the external 
gauge field (Aext), and the induced internal gauge field 
■A-int- Within single holon condensation scenario and us- 
ing loffe-Larkin formula, it can be shown that when both 
spinous and holons condense, we have Aint ~ —xAext- 
Since spinous are electrically neutral, they only couple to 
the internal gauge field and therefore they see —xA^xti 
so their efl^ective electric charge is —xe. Now we can 
estimate the value of a by computing the Green's func- 
tion of real electrons. Green's function can be calculated 
by convoluting the Green's function of holons with that 
of spinons. Within the single holon condensation sce- 
nario ge{k + eAex/c,uj) = Xcgs{k + eAint/c,u), where 
Xc is the condensation fraction of holon gas. This ex- 
pression means that the physical quasiparticles only re- 
sponse to the induced internal field Aint ^ —xA^xt, and 
therefore their effective electric charge is —xe. Therefore 
they carry xev^ current and this leads to a = x. On 
the other hand, in the double condensation scenario, as 
it has been discussed by Lee and Wen [21I, , we have 
ge{k + eAext/c,i^) = XpcQs {k + eAext/c,u}). Therefore, 
in this case, quasiparticles see the whole external field 
and we obtain a = 1 in consistent with the linear tem- 
perature coefficient of superfluid density measurements. 
Let us do more serious calculations now. In the pair 
condensed scenario, since Xpc (T) holons per lattice site 
condense at the ground state with energy Eg (which we 
showed before is exponentially small), we can write the 
diagonal part of the holons Green's function as 



igh [k, uj) = XpcS {k) 6{uj - Eg) + igh (k, lo) (27) 

where g^ denotes the uncondensed part of the system. 
The Green's function of spinons is 



igs {k,uj) = 
iAs (fc, Lo) — 



\Vk\ 



UJ + Es^k — «0+ UJ - Es^k — iO^ 

UkVk UkVk 



UJ + Es.k -iO+ UJ- Es.k - ^0^ 



(28) 
(29) 



Now let us choose the Coulomb gauge in which Aq = 
and V.A = 0. Since holons are charged quasiparticles, 
they couple to the both internal and external gauge fields. 
Since we have assumed pair condensation, the diagonal 



part of Green's function responds to the gauge fields and 
we have 

igA,h (fc, uj) = gh{k - e {Aint + A^xt) /c, i:^) (30) 

In the presence of the gauge field only the diagonal part 
of the spinons Green's function responds to the gauge 
field and the off-diagonal part does not change. Since 
spinons are electrically neutral, they only couple to the 
internal gauge field 

gA,s{k,uj) = gs{k - eAint/c,uj) (31) 

From = f} „hi it can be read that the Green's 

function of the real electrons is related to the Green's 
function of holons and spinons by convolution. 



gA,e {k,uj) = i 



gA,s{k + Q,uj + n)gA,h{Q,i^) (32) 



Let us separate the coherent and incoherent parts of 
the Green's function ge {k,uj) = 3^°'' {k,uj) + gl"'^ {k,uj). 
For the coherent part we immediately conclude 
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'!l{k,uj) = gl°^ {k + eA,xtlc,uo) 



(33) 



Which clearly implies that the effective charge of quasi- 
particles is — e and therefore they carry ev^ current, not 
xeVp . So we conclude that in the pair condensation sce- 
nario, quasiparticles carry the whole current and there- 
fore a = 1 in this case. 

In the single condensation scenario, the diagonal part 
of the holons Green's function is also igh{k,uj) = 
\{h)\ 5 {k) 5{uj)+ igh {k,uj). However the first term does 
not respond to the gauge field, like the offdiagonal part of 
the electrons Green's function and we have igA,h ik,uj) — 
\{h)fS{k)5{uj) + igh {k-e {Aint +Aext)/c,uj) . After 
convolution it can be verified that the coherent part 
of the real electrons Green's function is (7™e (^i ^) ~ 
gcoh j-^ _ QAint/c,uj). In the previous section we showed 
that in the SG state and at very low temperatures, Aint — 
-xAext, so we have g'j^^ {k,uj) = gl°^ {k + xeAext/ c,uj). 
This results a = x and therefore quasiparticles carry only 
X fraction of the eVp current. So the single holon conden- 
sation scenario gives x^ dependence for the linear tem- 
perature dependence of the superfluid density, which is 
far from experimental observations. 
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